^ 



Oh' 
< 



DESY 06-036 |hep-ph/060'4053 

DCPT/06/40, IPPP/06/20 
April 2006 



Next-to-Next-to-Leading Order Evolution 



o : of Non-Singlet Fragmentation Functions 



A. Mitov'', S. Moch'^ and A. Vogt^ 



> 

m 

O ' "Deutsches Elektronensynchrotron DESY 

(^ ' Platanenallee 6, D-15735 Zeuthen, Germany 

- "IPPP, Department of Physics, University of Durham 

""q,' South Road, Durham DHl 3LE, United Kingdom 

> 

H ' Abstract 

c^ ! 

We have investigated the next-to-next-to-leading order (NNLO) corrections to inclusive hadron 
production in e+e annihilation and the related parton fragmentation distributions, the 'time-like' 
counterparts of the 'space-like' deep-inelastic structure functions and parton densities. We have 
re-derived the corresponding second-order coefficient functions in massless perturbative QCD, 
which so far had been calculated only by one group. Moreover we present, for the first time, the 
third-order splitting functions governing the NNLO evolution of flavour non-singlet fragmentation 
distributions. These results have been obtained by two independent methods relating time-like 
quantities to calculations performed in deep-inelastic scattering. We briefly illustrate the numerical 
size of the NNLO corrections, and make a prediction for the difference of the yet unknown time- 
like and space-like splitting functions at the fourth order in the strong coupling constant. 



In this letter we address the evolution of the parton fragmentation distributions D and the 
corresponding fragmentation functions F^ in e+e^ annihilation, e+e^ — > y, Z ^ /i + X where 
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, , ^ „,l+cos^e)Ff + -sin^0F/^ + -cos0Ff . (1) 

Here represents the angle (in the center-of-mass frame) between the incoming electron beam and 
the hadron h observed with four-momentum p, and the scaling variable reads x = 2pq/Q^ where 
q with ^^ = 2^ > is the momentum of the virtual gauge boson. The transverse (T), longitudinal 
(L) and asymmetric (A) fragmentation functions in Eq. ([T]) have been measured especially at LEP, 
see Ref. [1] for a general overview. Disregarding corrections suppressed by inverse powers of Q^, 
these observables are related to the universal fragmentation distributions D^ by 

Fa\^,Q')= I f'^-^c^,{zMQ'))DH-.Q'). (2) 

The coefficient functions c^ f in Eq. @ have been calculated by Rijken and van Neerven in Refs. 
[2-4] up to the next-to-next-to-leading order (NNLO) for Eq. ([T]), i.e., the second order in the 
strong coupling a^ = as(2^)/(47t). Below we will present the results of a re-calculation of these 
functions by two approaches differing from that employed in Refs. [2-A]. 

Besides the second-order coefficient functions, a complete NNLO description also requires 
the third-order contributions to the splitting functions (so far calculated only up to the second 
order [5-7]) governing the scale dependence (evolution) of the parton fragmentation distributions. 
In a notation covering both the (time-like q, a = 1) fragmentation distributions and the (space- 
like q, Q^ = —q^, a = — 1) parton distributions, the flavour non-singlet evolution equations read 



d 



dlnQ^ 
with 



/-(x,e2) = f' ^p- (,,a,(e2)) /-(^,e2) o) 

Jx Z ^Z ' 



P-(x,as(e2)) = asP(0)"^(x) +a2/'ii^"^(x) +a3/'(2)-(^) + ... . (4) 

The superscript 'ns' in Eqs. Q and Q stands for any of the following three types of combinations 
of (parton or fragmentation) quark distributions, 

fik = li^^i - i<lk±qk) , r = Lr=lilr-qr) , (5) 

where rir denotes the number of active (effectively massless) flavours. As detailed below, we have 
obtained the so far unknown time-like NNLO splitting functions i'^J"^(jc) in Eq. ©. 

As already indicated in Eq. ©, the space-like and time-like non-singlet splitting functions are 
identical at the leading order (LO) [8], a fact known as the Gribov-Lipatov relation. This rela- 
tion does not hold beyond LO in the usual MS scheme adopted also in this letter. However, the 
space-like and time-like cases are related by an analytic continuation in x, as shown in detailed dia- 
grammatic analyses [5,9] at order a^, see also Refs. [10, 11]. Moreover, another approach relating 
the non-singlet splitting functions has been proposed in Ref. [12]. Hence it should be possible to 
derive time-like quantities from the space-like results computed to order a^^ in Refs. [13-15]. 
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We start the analytic continuation from the unrenormalized (and unfactorized) partonic trans- 
verse structure function F^ in deep-inelastic scattering, y*^ — > X (and correspondingly for Fi and 
^3 -^ Fa), calculated in dimensional regularization with D = 4 — 2e and the scale fj [13-15], 



F,\alQ^) = 5(1 -x) + l^ia^y ( ^ ) F^^, . 
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Q' 



-/e 



(6) 



The bare and renormalized coupling a^ and ttj are related by (recall a^ = as/(47r)) 
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0^2 



P^ Pi 



as <+ I ^-— |< -F ... 



(7) 



with Po = 1 1/3 C^ — 2/3 rir etc. The expansion coefficients in Eq. ^ are then decomposed into 
form-factor (^i) and real-emission (^/, defined analogous to Eq. ©) contributions [16] 
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2^2^(1 -x) + {!Fif8{l-x) + 2!Fi!lli+ 111 
2iF36(l-x)+2iFiiF25(l-x) + (2j2 + (i?'i)^)5?.i+2yi^2 + ^3 



(8) 



The analytic continuation of the form factor to the time-like case is known. The x-dependent 
functions Hii are continued from x to 1 /jc [5 , 9] , taking into account the (complex) continuation of q^ 
(see Eq. (4.1) of Ref. [16]) and the additional prefactor x^^^^ originating from the phase space of 
the detected parton in the time-like case [3]. Practically this continuation has been performed using 
routines for harmonic polylogarithms [17, 18] implemented in FORM [19]. The only subtle point 
in the analytic continuations is the treatment of logarithmic singularities for x ^ 1, cf. Ref. [9], 
starting with 

ln(l-;c) -^ ln{l-x)-\nx + i%. (9) 

Finally the bare transverse fragmentation function Fj is re-assembled analogous to Eq. (jSJ), 
keeping the real parts of the continued !^/ only, and the time-like non-singlet splitting functions 
and coefficient functions can be read off iteratively from the non- singlet mass factorization formula 
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2P^^l,+3P^l,c^^ +/'(0)(64') -3P(°)41) -3Po4^ 
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+ ... 



where the expansion coefficients Fj ^ now refer to an expansion in the renormalized coupling at the 
scale Q^ = q^. The products of the x-dependent (generalized) functions in Eq. (ITOI) are to be read 
as Mellin convolutions or, more conveniently, as products in Mellin-A^ space, employing routines 
for harmonic sums and their inverse Mellin transform back to .^-space [18-20]. 



Unlike the diagrammatic treatment of Refs. [5,9] — which cannot be emulated at order a^^ 
using the only available space-like calculation based on the optical theorem [13-15] — the above 
procedure is not entirely rigorous, as Eq. (HJ) does not represent a full decomposition according to 
the number of emitted partons. This would be required, for instance, in a subtraction formalism 
for exclusive observables at NNLO [21] and beyond. In Eq. © the real functions ^„>2 do not 
only collect, e.g., n-gluon tree-level amplitudes, but also combinations of real emission and virtual 
corrections. Especially, starting at order a^^, the decomposition Eq. ^ includes overlapping diver- 
gences from triple unresolved configurations when two particles become soft and one is coUinear. 
Thus one has to be prepared for some problem in the abelian (Cp) piece related to n^ contributions 
originating from phase space integrations over unresolved regions. 

At the second order, however, the above procedure works perfectly at least for the terms in the 
£-expansion written down in the second line of Eq. (ITOb . thus including the previously unknown 

(2) 

e-coefflcients aa {x), a = T, L,A, which we will not write out here for brevity. We have verified 
this by comparing to a direct calculation, to be presented elsewhere, of e^e^ ^ y, Z ^ h + X 
to this accuracy in e using the approach of Ref. [22]. Especially, we also re-derive the 0{a^) 

(2) 

coefficient functions Ca {x) which so far were only calculated in Refs. [2-4]. The differences 
between the time-like non-singlet coefficient functions for Eq. O and the corresponding (by the 
structure of the respective hadronic tensors) quantities in deep-inelastic scattering read 



(2) ^.) J2) 
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'^T.m^^ 



^Lnsi^) 



+ Q (^72- 8^2 + 176/3Ho-48Ho,o + 8H2-72Hi,o + (l+x) [-493/6+ 12^3 
+ IOC2 - 155/6 Ho + I6H0C2 - 6Ho,o + 12Ho,o,o - 2OH3 - IOH2 - 4H2,o - 29Hi 



-f42Hin-12H,i 



+ PqqW -84^3 -106/3^2 + 389/6 Ho -44H0C2 



+ 196/3 Ho,o - 24Ho,o,o + 36H3 + 62H2 - 20H2,o + 48H2,i - 40Hi ^2 + 134/3 Hi,o 
+ 32Hi,o,o + 56Hi,2 + 40Hij,o]+ 5(1 -x) [608/3^2 -24^2 

+ C^(Q-2C^)(^-76/3xHo + (l+;c)[-215/6 + 49/6Ho-9HiJ+pqq(x)[l2^3 
- 44/3 ^2 + 445/6 Ho - 2OH0C2 + 44/3 Ho,o + 24Ho,o,o + 4H3 + 22H2 - 4H2,o 



16Hi^2 + 22/3 Hi,o + 8Hi,2 - 8H1 j,o 



24H 



0,0,0 



+ 5(1 -jc) 466/3^2-24^. 



+ Pqq{-x) 32H_i,o,o + 8HoC2+16H_2,o 



+ CpnrU/3xlio + {l+x) 19/3 + l/3Ho + 2Hi +Pqq(x) 8/3^2-35/3Ho 



8/3Ho,o-4H2-4/3Hi,o +5(1 -x) -76/3^2 



(11) 



.(2) 



'^L,nsWU+e- 



•-L.nsWIep 



+ C^('583/9-860/9x- (14/3 + 152/3x)Ho- (12 +16x)Ho,o + (74/3 -112/3x)Hi 

-16(l+x)Hi,o + 8(l-2x)Hij+4(l+6x)[^2-H2]) 
+ C^(Q-2C^)('2317/45 + 8/5;c"i -3752/45^ + 32/5x2 + 32x^3 + (40 + 8/5;c"2)H_io 



- 1/5(98/3 + 8x^1 +728/3x-32x2)Ho + 16Ho-Lo + 32(x- 1/5x3) [^2 + H_Lo-Ho,o 
-8(1+2x)[h_i^2 + 2H_i,_i,o-H_i,o,o]+(1-2x)[46/3Hi-8HiC2 + 8Hi^o,o]) 
+ C^«^(^- 74/9 +112/9x + 4/3(l+4x)Ho- 4/3 (l-2x)Hi) , (12) 

+ C^((l-x)[-16C2 + 46Ho + 20Ho,o+16H2 + 24Hi,o]) + C^n;(^-4(1 -x)Ho 
+ Cp (Q - 2Cp) (;5qq(-x) [ - 64H_i,o,o - I6H0C2 - 32H_2,o + 48Ho,o,o 
+ 54(l-x)Ho + 32(l+jc)Hon 



.(2) 



(13) 

Here our notation for the harmonic polylogarithms //,ni,...,m»„(^)> "^; = 0, ±1 follows Ref. [17]. 
Furthermore we have employed the short-hand notation 



^o,...,o,±i,o,...,o,±i,...W 



H 



±(m+l),±(«+l), 



(14) 



suppressed the argument of the polylogarithms, and used the function Pqq{x) =2(1— x) ^ — l—x. 
The divergences forx — ;> 1 in Eq. dTTTl are to be read as plus-distributions. 

Eqs. (fTTT) - (fOl agree with the results in Refs. [3,4] up to a few (presumably typographical) er- 
rors in those articles. Specifically, forcj^ (z) the term —3Cp (1 +z)ln^zinEq. (A.6) of Ref. [3] has 
to be replaced by -3C^ (l+z)ln3z, and the contribution 2A{C}-Cjfp/2) \nz/{5z^) inEq. (A.8) 
by 24 {Cp — C^Cp/2) lnz/(5z) . The argument of Si,2 should read —z instead of 1 — z in the first 
term of Eq. (A.15) for cf ^ (z). Finally, in Eq. (17) of Ref. [4] for cf (z) the term 24C2 lnz/(5z2) 
has to be replaced by 24C^ lnz/(5z). We have also re-calculated the second-order gluon and 
pure-singlet coefficient functions for Cj and c^, finding complete agreement with Ref. [4]. 

We now turn to the corresponding NNLO (third-order) splitting functions Z'^J"^ in Eq. (HI). 
Using the same notation as above these functions are given by 



e\"w 



d(2)+ , 

CT=-P 



5/'(2)+(;c) 

+ \6CJ (pqq [x] [3 1 1 /24 Ho + 4/3 H0C2 - 1 69/9 Ho,o + 8Ho,oC2 - 22Ho,o,o 

- 268/9 Hi,o + 8Hi,oC2 - 44/3 Hi,o,o - 268/9 H2 + 8H2C2 - 44/3 H2,o - 44/3 H3 
+ (l+x)[-4Ho,oC2 + 25/2Ho,o,o + H2,o + 2H3]-(l-x)[325/18Ho + 50/3Hi,o 
+ 50/3 H2 + (3 -5jc)HoC2 - (173/18 - 691/18x)Ho,o) 

+ 16C^(Q-2C^)(pqq(x)[l51/24Ho + HoC3 + 13/6HoC2-169/18Ho,o + 8Ho,oC2 

- 13/2 Ho,o,o-8Ho,o,o,o- 134/9 Hi^o + 4Hi^oC2- 22/3 Hi,o,o-6Hi,oAO- 134/9 H2 



8H 



-3,0 



+ 4H2C2 - 22/3 H2,o - 2H2,o,o - 22/3 H3 - 2H3,o - 6H4J + Pqq(-x) 

+ 8H_2C2 + 8H_2,-i,o + 3H_2,o - 14H_2,o,o - 4H_2,2 + 8H_i,_2,o + 16H_i,_i,o,o 

+ 8H_i,oC2 + 6H_i,o,o - 18H_i,o,o,o - 4H_i,2,o - 8H_i,3 - 7H0C3 + 3/2HoC2 



- 8Ho,oC2 - 9/2 Ho,o,o + 8Ho,o,o,o + 2H3,o + 6H4J - ( 1 + x) [4H_2,o + 8H_ 1 ,0,0 
+ {l-x) [4H_3,o + 4H_2,o,o - 88/9 Ho + SHqCs - 28/3 Hi,o - 28/3 H2] - AxUq^i 

- (50/9 -184/9x)Ho,o-4;cHo,oC2 + (l 1/2 + 35/2x)Ho,o,o + 8xHo,o,o,o) 
+ 16C^n^(pqqW[-ll/12Ho-2/3HoC2 + ll/9Ho,o + 2Ho,o,o + 20/9Hi,o 

+ 4/3Hi,o,o + 20/9H2 + 4/3H2,o + 4/3H3]-(l+x)Ho,o,o + (l-^)[l3/9Ho 

+ 4/3Hi,o + 4/3H2]+(8/9-28/9x)Ho,o) , (15) 

and 

+ 16C| (Q - 2Cp) (;?qq(-x) [l6H_3,o - 16H_2C2 - 16H_2,-i,o - 6H_2,o + 28H_2A0 
+ 8H_2,2 - 16H_i,_2,o - 32H_i,_i,o,o - 16H_i,oC2 - 12H_i,o,o + 36H_i,oao + 8H_i,2,o 
+ 16H_i,3 + I4H0C3 - 3H0C2 + 16Ho,oC2 + 9Ho,o,o - 16Ho,o,o,o - 4H3,o - I2H4 



+ {\+x) 8H_2,o + 16H^i,o,o + 8HoC2-4H2,o-8H3 -{l-x 



8H_3,o + 8H_2AO + 10Ho 



+ 6H0C3 + 4Ho,oC2 - 8Ho,o,o,o + 8Hi,o + 8H2 - ( 10 - 6x)Ho,o - (12 + 24x)Ho,o,o 



(16) 



for both ^ = — and ^ = v. Eq. (15) is, in fact, not quite the result of the analytic continuation as 
described above, which returns a different coefficient for the term Cp pqq (x) Ho,o C2 in the first line. 

We have corrected this term by imposing the correct (vanishing) first moment of /'^Jj • As it is 
conceivable that the present form of the analytic continuation leads to other problems not affecting 
the first moment — recall the discussion in the paragraph below Eq. (ITOb — we obviously need a 
second, independent confirmation of our new results (ITSt and (IT6b . 

For this purpose we adopt the approach of Dokshitzer, Marchesini and Salam [12] (see also 
Appendix Bl of Ref. [23] ), where Eq. Q is rewritten as^ 



dlnQ- 



J X Z ^Z 



(17) 



and the modified splitting functions P^^^ are postulated to be identical for the time-like and space- 
like cases. Working out the perturbative expansion of the integrand along the lines of Ref. [12] one 
arrives at a successful 'postdiction' for the NLO difference P^Ji^{x) — P^J^\{x) of Refs. [5,7] 
(see Eq. (4) of Ref. [12] ) and the new NNLO prediction (^ = +, -, v, recall Eq. ©) 



Inx-P^i)^ 



®P^^^ + 



with ® denoting the Mellin convolution (cf. Eqs. Q and Q) and 



Injc-P^'^)] ^/'(i)^ 



(18) 



(19) 



'Ref. [12] also includes a shift in the argument of as in Eq. ( I17> which is irrelevant for our purpose of relating the 
time-like and space-like results. Note also that the notation for the as expansion in Ref. [12] differs from Eq. (0}. 



The evaluation of Eq. (ITSb . again performing the convolutions via a transformation to Mellin-A^ 
space, yields exactly Eqs. (fTSt and (fT6t . thus providing both the desired confirmation of these 
results and further evidence supporting the ansatz (flTt . 

Consequently, it is possible to make even a prediction for the fourth-order (N^LO) difference 
§p(3)^ of the (both unknown) time-like and space-like non-singlet splitting functions on this basis. 
Using the notation of Eq. (fTSt together with A®^ =Ai^A etc, this prediction reads 

5i'(3)^(^) = 2|[ln;c-p(2)^]^p(0)^[ln;c-/'(°)]®p(2)^^[lnx-p(i)^]®p(i)^} 

-2/'(0)®[ln;c-/'(0)]®3_4[/'(0)]®2®[lnjc-/'(«)]®[ln2x-/'(0)] 
_2/3[p(0)]®3^[^j^3_^.p(0)] _ (20) 

For brevity we refrain from writing out the resulting explicit expressions which are, of course, 
much more lengthy than Eq. (ITSl) . We expect that a first check of Eq. (l20b . or rather its first 
line (which dominates the large-x behaviour, cf. Ref. [12]), will be obtained via next-to-leading 
order calculations in the large-n . expansion, generalizing the leading-n ^ result of Ref. [24] . The 
corresponding contribution to Eq. (OUt is identical for all three non-singlet cases and reads 

16/81C^«/(77qqW[- (159/4- 120^2 + 36C3)Ho + (23 +72C2)Ho,o-279Ho,o,o 
- 216Ho,o,o,o - 76(Hi,o + H2) - 240(Hi,o,o + H2,o + H3) - 108(Hi,o,o,o + H2,o,o 
+ H3,o + H4)] + (l-x)[-(260-72C2)Ho-276(Hi,o + H2)-144(Hi,o,o + H2,o 
+ H3)] + (l+x)108Ho,o,o,o-(466-398x)Ho,o-(90-450.^)Ho,o,o) • (21) 



The a"Cpn'r contributions dominating bP^" '^'"** in the large-«^ limit are given to all higher 
orders n by a straightforward generalization of Eq. (ITHl and the first line of Eq. (EUt . 

Returning to the NNLO coefficient-function and splitting-function differences dTTT) - (IT6b . we 
note that these functions include harmonic polylogarithms up to the same weights as the corre- 
sponding space-like results [13-15,18], with the interesting exception that weight-4 functions enter 
Eq. (1131 only for the SU(nc) group-factor combination Ca — 2Cf suppressed as l/ric in the limit of 
a large number of colours Hc . Except for the longitudinal coefficient function CL,ns , the differences 
between the time-like and space-like quantities are parametrically suppressed in the large-jc limit, 
since the logarithmically enhanced soft-emission contributions to 5?,/ in Eq. (jSJ) are, as they have to 
be, invariant under the analytic continuation. Specifically, the (identical) leading large-jc contribu- 
tions for Cj ns — Cj j(g and c^ ns ~''3 ns contain plus-distributions only up to [(1 — ;c)^Un(l —x)] + 
(all proportional to TT^-terms arising from the analytic continuation of the form factor), and those 
for the splitting-function differences read, 

6p(2),ns^^^ J) = -4Aj^^^hn(l-x) + 0(1) , (22) 

as predicted in Ref. [12], where Aq are the coefficients of a" [1 —x]^ in Eq. Q, cf. Ref. [13]. 
The leading small-x terms, on the other hand, differ between the space-like and time-like cases. 



The numerical impact of the second-order contributions to the time-like coefficient functions 

(2) 

Ca {x),a = r, L, A, has been discussed in some detail already in Refs. [2-4]. Here we confine our- 
selves to the transverse fragmentation function Ft, the largest contribution to the right-hand- side 
of Eq. ([T]), see Ref. [1]. In Fig. [T] we compare the corresponding non-singlet coefficient function 
Cj^^ to its counterpart Cj ^^ in deep-inelastic scattering. In order to facilitate a direct comparison, 
the same schematic shape has been used for the non-singlet fragmentation distributions and parton 
distributions. As obvious from the figure, the higher-order corrections for Cj are in general con- 
siderably larger than those for Cj. Nevertheless the second-order term changes the NLO results, 
under the conditions of Fig.[Tl by 5% or less from x ~ 0.55 down to very small values of ;c. 

The pattern is quite different for the non-singlet splitting functions illustrated, in a similar 
manner but at a lower scale, in Fig. El In A'^-space, for example, the ratio 5/'"Y(/'^!^_j ~^lc)) 
quickly decreases with increasing N, at NNLO from about 1/2 at A^ = 2 to about 1/6 at A'' = 8. 
Consequently the total time-like splitting functions i'^ij is only mildly enhanced, e.g., by 8% and 
2% for these two values of N and tts = 0.2, with respect to their space-like counterparts ^'^=-1 
discussed in detail in Ref. [13]. As shown in the right part of Fig.|2l the small-x scaling violations 
of the non-singlet fragmentation distributions are weaker than those of the parton distributions. 
For the chosen input distribution this reduction increases, in a perturbatively stable manner, from 
about 10% at X = 10"^ ^^ ^bout 30% at x = 10"^. 

To summarize, we have re-derived the O(cXs) coefficient functions [2-4] for the inclusive 
production of single hadrons in e+e annihilation [1] and obtained, for the first time, the cor- 
responding third-order splitting functions for the flavour non-singlet fragmentation distributions. 
Our derivation of the latter quantities rests on relations between the time-like and space-like cases, 
see especially Refs. [5,7,9, 12], and the third-order calculation of deep-inelastic scattering of 
Ref. [13-15]. We expect that a further study of these relations, backed up by fixed Mellin-A'^ cal- 
culations along the lines of Ref. [22], will facilitate an extension of our derivation to the NNLO 
flavour-singlet splitting functions and, at least for Fl, the (a^ ) coefficient functions. 

Once this step has been taken, the way is open for full NNLO analyses, e.g., along the lines 
of Ref. [25], of high-precision data on e^e~ ^ h + X from LEP and a future International Lin- 
ear Collider. Due to the universality of the splitting functions, our results also represents a first 
step towards NNLO analyses of high-pr hadron production in ep and pp collisions, where very 
large NLO corrections strongly suggest sizeable higher-order contributions, see Refs. [26-28] and 
[29-31], respectively. Another application concerns the Z7-quark spectrum in top decays (see 
Ref. [32]) where, after the calculations of Refs. [33,34], the third-order time-like splitting func- 
tions will facilitate a complete NNLO treatment in the framework of perturbative fragmentation. 

Form and Fortran files of our results can be obtained from |http: //arXiv.orql by down- 
loading the source of this article. Furthermore they are available from the authors upon request. 

Acknowledgments: We are grateful to L. Dixon, R Uwer and W. Vogelsang for useful discus- 
sions. Our numerical results have been computed using the Fortran package of Ref. [35]. A.M. 
acknowledges support by the Alexander von Humboldt Foundation. The work of S.M. has been 
supported in part by the Helmholtz Gemeinschaft under contract VH-NG-105. 
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Figure 1 : Comparison of the coefficient function Cj ^^ for the (time-like) process e+e ^ /i + X 
with its counterpart Cj j^^ in (space-like) deep-inelastic scattering for Q^ ~ M\ . Left plot: Mellin 
moments, right plot: convolutions ^ with a schematic input shape denoted by /. 
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Figure 2: The differences dP^^ = P^^i — P^^_i between the time-like (o = 1) and space-like 
(o = —1) non-singlet splitting functions at a 'low' scale characterized by the (order-independent) 
value tts = 0.2 of the strong coupling constant. Left: moment-space comparison with the higher- 
order corrections in the space-like case. Right: convolutions with a schematic input shape. 
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